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ABSTRACT 


The quest for the highest resolution microwaves imaging and the principle of 
time-domain imaging has been the primary motivation for recent developments in 
time-domain techniques. With the present new technology fast time varying 
signals can now be measured and recorded both in magnitude and in phase. It has 
also enhanced our ability to extract relevant details concerning the scattering 
object. 

In the past, the inference of object geometry or shape from scattered 
signals has received substantial attention in radar technology. Various inverse 
scattering theories were proposed to develop analytical solutions to this 
problem. Furthermore, the random inversion, frequency swept holography, and the 
synthetic radar imaging, all of which have two things in common: (a) the 
physical optic far- field approximation and (b) the utilization of the channels as 
an extra physical dimension, were also advanced significantly. 

Despite the inherent vectorial nature of electromegnetic waves, these scalar 
treatments have brought forth some promising results in practice with notable 
examples in subsurface and structure sounding. The use of time-domain imaging 
for space robotic vision applications has been proposed. A multi sensor approacah 
to vision has been shown to have several advantages over video-only approach. 
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OBJECTIVES 

The continued research efforts for this project should be as follows: (1) 

develop the application of time-domain imaging for robotic vision, in order to 
emphasize this we stress that it is the scattered field not the total power that 
is of interest; (2) research the hardware/technology needed for this application; 
(3) derive theoretical expressions for planar, cylindrical, spherical, and 
octobox configurations (the images will be synthesized on the basis of 
time-domain reflectance); (4) study the effects of pulse width, polarization, 
look angle, and carrier frequency and (5) compute magnitude and phase images. 
The limitations of the present analytical techniques will be identified. Also, 
experimental verification of the analytical results will be performed to develop 
hardware approaches. 


INTRODUCTION 

The detailed literature review conducted this summer at the Johnson Space 
Center, in the area of Time-domain microwave target imaging, led to a research 
topic of determining the size, shape, and electromagnetic properties of a 
scatter, given the incident and scattered eletromagnetic fields. The problem of 
interest here is the reconstruction of the shape of a convex scatter from 
knowledge of high frequency far field scattered from the object in response to a 
known incident field. 

The quest for the highest resolution microwaves imaging and the principle of 
time-domain imaging has been the primary motivation for recent developments in 
time-domain techniques. In the last decade, modern techniques in sampling 
devices and advent of fast pulse generators has brought new technology to the 
practicing radar engineers in the measurements of picosecond at a much reduced 
cost. With the present new technology fast time varying signals can now be 
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measured and recorded both in magnitude and in phase. 

Todays technology has enhanced our ability to extract relevant details 
concerning the scattering object, which conventional radar ranging lacks. 

In the past, the inference of object geometry or shape from scattered 
signals has received substantial attention in radar technology largely due to its 
adademic significance and the understandable commercial and military values. In 
the last decade various inverse scattering theories were proposed to develop 
analytical solutions to this problem. Furthermore, the random inversion, 
frequency swept holography, and the synthetic radar imaging, all of which have 
two things in common: 

(a) The physical optic far-field approximation. 

(b) The utilization of the channels as an extra physical dimension, were 
also advanced significantly. 

Despite the inheret vectorial nature of electromagnetic waves, these scalar 
treatments have brought forth some promising results in practice with notable 
examples in subsurface and structure sounding. Recently Dr. Kumar Krishen.of 
NASA/JSC, has proposed the use of time-domain imaging for space robotic vision 
applications. A multisensor approach to vision has been shown to have several 
advantages over video-only approach. 


PROBLEM 

(THEORETICAL APPROACH) 

The problem of scattering from rough surfaces is inherently different in 
nature from that of scattering by other bodies. Usually in the rough-surface 
problem, an exact knowledge of the shape of the surface is neither available, nor 
is it of interest to the radar operator. Instead, only average properties of the 
surface shape enter into the problem. It is clearly understood that the last 
requirement rule out a boundry-val ue approach, since the exact boundry is not 
known. We are more interested in the relationship between the average scattered 
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field or radar cross section and the average surface properties. The radar cross 
section is normalized by the area A, and defined as, 
7 (S t ,e s ,9, ) = <cr (c ^ ,c s ,p s > /A, where the three angles 

S. ,6 , and o are shown for bistatic case in figure la below. The rough 

US S 

surface is assumed to consist of height variations about a mean plane, which is 
taken as the xy-plane. The incident wave lies in the xz-plane at a polar angle.. 


Z 



Figure la - BISTATIC SCATTERING GEOMETRY 



Figure lb - THE "SCATTERING GEOMETRY" 

Figure la shows the bi static scattering for a planar rough surface whose mean 
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height coincides with the xy-plane and figure lb shows the scattering geometry. 
I is the plane of incidence and S is the scattering plane. 

The < > refer to an average. For the normaized average backscattering 

cross section are defined as 


) - y (S ,6 

Generally one has specific purposes or reasons for investigating scattering 
from a rough surface broadly these ultimate goals may be divided into three 
categories: 

(1) The problem of direct scattering; one wishes to know the average 
properties of the scattered signal or cross section when the surface properties 
of the rough surface are known, and the scattering information is expressed in 
terms of the surface properties. 

(2) The problem of inverse scattering; one wishes to obtain statistical 
information about the rough surfaces from a knowledge of the average properties 
of the scattered field. This problem is more difficult in that there appear to 
be many classes of rough surfaces producing the same average scattering cross 
section as a function of the bistatic scattering angle and wavelength. 

The (scattered field) radar cross section, can be expressed as nr times the 
power delivered per unit solid angle in the direction of the receiver divided by 
the power per unit area incident at the target. The factor nr enters from the 
definition of solid angle. Assuming for the moment that the propagation path 
between the target and the receiving system is lossless, then this power ratio 
may be expressed as 

4t : ( 2* ♦ 2* # ) / ( Z 1, ♦2*' * ) =* 4-r : ( H* ♦ E* * )/(H l ) 
where E s and H s are the scattered electric and magnetic fields, respectively and 
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E" and h*' are the incident fields. The scattered field is defined to be the 
difference between the total field (with the targety present) and the incident 
field E 1 (the total field that would exist if there were no target present). 
This is summarized as E S * E - E 1 . 

Finding the cross section (a) now becomes a problem in electromagnetic 
field theory. In order that the cross section be independent of r, it is 
desirable to let r in (Eg., 3) become arbitrarily large (scalar definition, not a 
tensor function) . 


= 4-^g r : (E J • Z* *)/{?*?* ) = 4 r s ( H* * E* * ) / ( K 1 ' * H v ) 


To compute the scattered field one sometimes computes the current induced on 

the target and then treats the target current distribution in terms of an 

equivalent aperature distribution. Antennas are often compared to an isotropic 

antenna, that is, an antenna which radiates uniformly in all directions. 

In the clutter problem, return from terrain is not wanted. However, in many 

cases the presentor clutter is unavoidably along with the desired signal, one can 

detect and analyze the desired signal significantly if more is known about the 

properties of the clutter or noise produced by terrain scattering. It is noted 

that surface information is generally known, and the properties of the scattered 

* * 

signal are sought related to (Eq. 3), At this point in time, a few more 
mathematical definitions may be of interest. 

The source signal is related to the Fourier transform, F(w), by: 


(a) 7(w) - 




dt. 


(b) f(t) = j^rf X F(w)e- jw *dt 
* 


Known by the Fourier transform pair as follows: 


(a) ?(w) * /°°f ( t ) e J w ‘* dt , (to) f(t) = i-^/°F(w)e’ j wi ' di 

- OO - 00 
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Capitals are used to denote transformed function unless otherwise noted or 
stated. The wave number k = w/c, where c is the speed of light in the medium. 
a and pare the direction cosines of the vector r projected onto the xy-plane 

(with cc * sin? cosr , [3 * sin? sinS ) . 

It can be assumed that the field distribution over the antenna apertures and 

object plane is space - time separable; the field strength of each radiating 
element observed at some point r is e(r,t) » e(r)f(t— £•) 
the field strength. Where f(t- §-) is the sine wave traveling at the speed 
of light. 

Consider the scatter in figure 2 to be placed in a far field, such that the 
assumption of a uniform impinging planewave over the object can be assured. The 
incident wave form i(t), on the object surface is obtained by the application of 
the far field Ki rchhoff-Huygens principle, namely 



Figure 2.- SCHEMATIC DIAGRAM OF THE TARGET IMAGING SYSTEM 

Let P be a point of observation and let R’ be the distance from P to a point 
;(x,y) on the surface S. The scattered field E at P is given by the Helmholtz 
integral as follows: 
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5(P> 


» ^ CHf? Vft)ds 

Let V be a volume bounded by a surface S and let A be any continuous vector 
function. Then, according to the divergency theorem, 

JJJ 7* Adv =££ A * nds 


Green’s First Theorem: 


/// v 7E,T7 V dv ■ ££ E f? ds " /// v 2 vav 

Green's Second Theorem: 

/// V (37‘V7 ! E)dv -£$<*£?- V £| )ds 


Note, the scattered field E(p) at the point of observation is given by the 
Helmholtz integral as follows: 


2 fp) 



e^££)ds 


(To get a better mathematical understanding of the application on how the 
scattered wave form is obtained, please refer to the author on the derivation of 
the scattered wave of interest in frequency domain). 


I(w) * «wF(w)f fe(x, y)exo c ‘ r o : axav 

o -x 

I (w) = - ° f E *wF(w)T(w)2. (to, r kf3. }exo c " ' o : 

* “ e • e v * * 


Where 2. (kci, ,kj3 4 ) is the transmitting-antenna anisortopic directional 

V * * 
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filter and r is the distance between the transmitter aperture and an object plane 
taken as reference. 

By stationary - phase evaluation and a second application of the 
Kirchhoff-Huygens principle, it can be shown, with little algebra that the 
received waveform, g(t), is 

00 

g(t) « const/ F(w)T(w)R(w)[H (kct ,kp )+ 
i 00 c c ^ c ' 

(jw) 2 2, (k*. y ,k0. )2 ? (kct. ,k(3 f }2 0 (kct Q ,k(3 0 )e‘ J:<? dw 

For the monostatic case commonly encountered (see figure 3 below), the 
transmitter and receiver coincide. In particular, if the phase center of the 
object is on boresi.ght which suggests that the target is closely tracted, eqr.16 

f 

can be considerably simplified to yield: 

. 0 ° 

g(t) = const / F(w)T(w)R(w) [2 n (lei .kj3, ) + 

- oo a « c 

( jw) : E v ( 0 , 0 ) E f (0,0)2 e (k* 0 ,k(3 e )e' Jwi dw 

The delineation of the object response is exercised as follows: First, the 

time response in the absence of any scattering objects is recorded. 



Figure 3 - IMAGING GEOMETRY 
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THE ROTATION TRANSFORMATION • 



This signal is subsequently substracted from other records with the scatter 
present, thus removing the antennas coupling filter. Given a sufficiently large 
free space, multiple path scattering between the object and surrounding obstacles 
can be gated out in time. By definition 

z 0 (kst 0 ,kJ3 0 ) =* f felx., y)exo c 1 x ‘ *0 * * ,!3 o y * : axdy 
- X 

Having removed the coupling term, the backscattering waveform g(t), gives the 
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system kernel, K(w), 

:<(w) = cost+(jw) : F(w)T(w)R<w) 2 . ( 0 , 0 ) 2 , ( 0 , 0 )e' : J;,? o 

with the use of mathematical techniques, the backscattering waveform was 
simplified, to give the system kernel, K(w), further- giving an opportunity to 
obtain the electric field strength in frequency domain 

and by the Projection theorem , the response obtained is equivalent to a slice in 
the two dimension fourier plane (Figure 4). 


I 



Figure 4. - REPRESENTATION OF THE ONE-DIMENSIONAL TARGET FREQUENCY RESPONSE 
IN THE TWO-DIMENSIONAL FOURIER PLANE AT A PARTICULAR VIEWING 
DIRECTION 


It is further noted that the help of figures 3a, 3b assisted in a change in 
spatial coordinates in Figure 3 from one dimension coordinate system to the two 
dimension coordinate system 2(ki,I<3) . 

Followed by an inverse fourier transform (Eq. 11) yields the object impulse 


response 


V " J *' **' "c* dw = e(t - 2x'sin|-> 


- oo 
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which can be considered as a series of scans across the object surface where each 
is an integral of the object field which is orthogonal to the viewing direction. 
The spatial resolution is thus achieved by the utilization of the time channel 
and the minimal resoluable dimension Ax is related to the sampling integral T, 
namely 


Ax = 


1 S V 


7F 


The two-dimensional object field distribution is recovered, namely the 
reconstruction field strength, e:(x,y) 

e(x, y) * -±-T f* /°Z(k£& f lC|3 ) A ( kct. , icJ3 )e : ‘ ftx ' |ly) kdxa? 

<rr 5 - oo 


a smoothing filter is introduced to compensate due to the factor k in 

the field strength equation e e (x,y) . It is not necessary to have A, for' 

noise free data. For practical numerical evaluation, the finite sampling rate 
will impose an upper limit, R, to the wave number k. To eliminate abrupt cutoff 
which would cause ringing in the resultant image, it is necessary to introduce 

A(k*,ld3) - siii(r^) / (|P) 

* • 

which yields a smoothed image 


e(x,v) - -l, 

4 r" o - 2 


and is the basis of our reconstruction field strength. 
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RESULTS 


Theoretical calculations of the single and train pulse of the schematic diagram 
(figure 2) of the target imaging sytem led me to analyze fourier transform of the 
input function. A program was written in Basic and the use of BOSS (Block 
Oriented System Simulator) enable me to plot the energy supplied to the target 
and give the time plot of a single and train of pulses. I was also able to 
verify the magnitude spectrum of a single pulse continuous descrete function and 
the magnitude spectrum of a line spectrum of pulse train. 


APPLICATIONS 

0 Recently NASA/Johnson Space Center has proposed the use of time-domain 

imaging for Space Robotic vision applications. 

0 Many space objects are covered with heat shields, etc. (Dielectric 

Materials) and microwave can sense through this layer. 

0 A multisensor approach to this vision has been shown to have several 
advantages over video-only approach. 

RECOMMENDATIONS 

0 It is recommended that the second phase of the research project, inititated 
at the NASA/Johnson space Center, in the area of Time-Domain Robotic Visions 
Applications, be continued to study the effects o9f pulse width, 

polarization, look angle, and phase images on robotic vision, and to 
conduct experiments to verify the theoretical results. 

o The amount of time provided for the theoretical and experimental 

verification of the Time-Domain Microwaive Robotic Vision Application was 
not sufficient in order to perform the above experimental verification. 
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LIST . 2 

5 PRINT "The Energy Supplied To The Target is lF( w )| = F ( w ) F ( w ) * " 

6 PRINT r 1 

7 PRINT 

8 PRINT 't 

9 PRINT 

10 FOR N = 1 TO 4 

15 PRINT “N" , "w" , "F(w)F(w)*" , " B", 

16 PRINT 

20 FOR W = 5 TO 10 

30 B = C 2/ ( N * 3.14 ))* (£t N (N * .0314)) * (W/( W*W - N*2*3 . 14 )/EZ?( 10-10 ) ) 

35FW=(l/(2* 3.14)* (.01))+ B * (1/(2*3.14) *(.01))-3 

60 PRINT N, W, FW, B 

70 NEXT W 

80 NEXT N 

O 


RUN 

The Energy Supplied To The Target 


is |F(w)|‘ 


= F(k)F(w)* 


|f(w)| = 


N 

w 

F(w)F(w)* 

B 

1 

5 

6 . 941864E-03 

5. 3410C3E-03 

1 

6 

5. 625806E-03 

4.027021E-03 

1 

7 

4 . 8741 69E-03 

3.27S615E-C3 

1 

8 

4 . 36E317E-03 

2.771547E-C3 

1 

9 

4.0047E9E-C3 

2 . 40859SE-C3 

1 

10 

3 . 729419E-03 

2 . 123665E-C3 

N 

w 

F(w)F(w)» 

B 

2 

5 

9.63845E-03 

8.033302E-03 

2 

6 

6 . 71S522E-03 

5 . 116069E-C3 

2 

7 

5.437B77E-Q3 

3 . S294C7E-C3 

2 

8 

4.7C5582E-03 

3.10837SE-03 

2 

9 

4 . 224854E-03 

2 , F 2E2: 2E-03 

2 

10 

3 . 881775E-C3 

2.28577SE-C3 

N 

w 

F(w)FCw)* 

B 

3 

5 

1.7B2794E-02 

1.E2Q377E-C2 

3 

6 

8 . 586145E-03 

6 . 982569E-C3 

3 

7 

6 . 234785E-C3 

4 . S2504SE-C3 

3 

8 

5.125711E-03 

3 . 537721E-C3 

3 

9 

4.488432E-03 

2.691472E-03 

3 

10 

.0040559 

2 . 4 60625E-03 

N 

w 

FCw)F(w)* 

B 

4 

5 

- . 8308692 

-.8311381 

4 

6 

1 . 261C31E-02 

1 . 100044E-C2 

4 

7 

7.448911E-03 

5 . 847244E-C3 

4 

8 

5 . 7C2307E-C3 

4.104415E-C3 

4 

9 

4.ei019SE-03 

3.21272SE-C3 

4 

10 

4 . 2SC522E-C3 

2. 663922E-C’3 


Table I - The Energy Supplied To the Target 
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Figure 8.- Magnitude Spectrum of A Pulse Train 
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EXPERIMENTAL VERIFICATION 



Figure 9.- Experimental Arrangement of the Target 
Scattering Measurement 

The above is a circuit diagram which will be used to verify the results of 
elliptical and bowtie geometry of targets (Figures 10, 11, 12, and 13) 
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Reconstructed Images -Ellipse 
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